/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

phases (particles gas);
particles
{
    phaseType        granular;
    diameterModel constant;
    constantDiameterCoeffs
    {
        d           10e-6;
    }

    alphaMax    0.63;
    alphaMinFriction 0.5;

    type basic;
    thermoType
    {
        equationOfState rhoConst;
        thermo eConst;
        transport constIso;
    }
    specie
    {
        molWeight  100.0;
    }
    equationOfState
    {
        rho0         2500;
    }
    transport
    {
        kappa       0.36;
    }
    thermodynamics
    {
        type        eConst;
        Cv          745;           // Heat capacity
        Hf          0.0;
    }

    kineticTheoryCoeffs
    {
        e           0.999;
    }
    residualAlpha 1e-10;
    residualRho 0;


}

gas
{
    phaseType       fluid;

    diameterModel constant;
    constantDiameterCoeffs
    {
        d           0;
    }

    type multicomponent;
    species      (He N2);
    defaultSpecie He;
    calculateDensity yes;
    thermoType
    {
        transport   const;
        thermo      eConst;
        equationOfState idealGas;
    }

    He
    {
        equationOfState
        {
            gamma           1.667;            // Heat capactiy ratio
        }
        specie
        {
            molWeight       4.0;
        }
        transport
        {
            mu              1.96e-5;
            Pr              1;
        }
        thermodynamics
        {
            Cv          3115.6;
            Hf          0;
        }
    }
    N2
    {
        equationOfState
        {
            gamma           1.4;            // Heat capactiy ratio
        }
        specie
        {
            molWeight       28.00;
        }
        transport
        {
            mu              1.76e-5;
            Pr              1;
        }
        thermodynamics
        {
            Cv          718;
            Hf          0;
        }
    }
    residualRho     1e-6;
    residualAlpha   1e-10;
}

blending
{
    default
    {
        type none;
        continuousPhase gas;
    }
}

interfacialPressure
(
    (gas and particles)
    {
        type single;
        phase gas;
    }
);

interfacialVelocity
(
    (gas and particles)
    {
        type single;
        phase particles;
    }
);

aspectRatio
(
);

solveDragODE off;
dragODECoeffs
{
    solver seulex;
    relTol 1e-3;
}

drag
(
    (particles in gas)
    {
        type GidaspowErgunWenYu;
        residualRe 1e-3;
        swarmCorrection
        {
            type none;
        }
    }
);

virtualMass
(
);

heatTransfer
(
    (particles in gas)
    {
        type RanzMarshall;
        NuModel Gunn;
    }
);

lift
(
);

wallLubrication
(
);

turbulentDispersion
(
);

kineticTheory
{
    alphaMax 0.63;
    residualAlpha 1e-8;
    alphaMinFriction 0.5;

    radialModel SinclairJackson;
    viscosityModel none;
    conductivityModel Syamlal;
    granularPressureModel Lun;
    frictionalStressModel JohnsonJackson;

    JohnsonJacksonCoeffs
    {
        alphaMinFriction        0.5;
        Fr                      0.05;
        eta                     2;
        p                       5;
        phi                     28.5;
        alphaDeltaMin           0.05;
    }
}
// ************************************************************************* //
